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In light of the proton radius puzzle, the discrepancy between measurements of the proton charge radius 
from muonic hydrogen and those from electronic hydrogen and electron-proton scattering measurements, we 
reexamine the charge radius extractions from electron scattering measurements. We provide a recommended 
value for the proton RMS charge radius, r® = 0.879 ± 0.011 fm, based on a global examination of elastic e-p 
scattering data. The uncertainties include contributions to account for tension between different data sets and 
inconsistencies between radii using different extraction procedures. 
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I. INTRODUCTION 

Multiple extractions of the proton RMS charge radius have been made based on global analyses of world data on elastic e-p 
scattering. The radius is related to the slope of the charge form factor, Ge(Q 2 ), at Q 2 = 0. As such, the extraction is sensitive to 
the precision of the data available at low Q 2 , the minimum Q 1 value of the data set, the uncertainty in the absolute normalization 
of the data, and the model dependence of the fitting function used to perform the extrapolation to Q 2 = 0. Many radius 
extractions have not examined these questions in detail or have not included uncertainties associated with the model dependence 
of the extraction, yielding a significant range of quoted values for the charge radius, often with unrealistic uncertainties. More 
recently, a high-precision data set from Mainz has significantly expanded the body of data at low Q 2 Oh- 

Extractions which attempt to account for all of these uncertainties appear in Refs. . All of these extractions yield proton 
RMS charge radius values of approximately 0.88 fm, consistent with those extracted from atomic physics measurements in 
electronic hydrogen @|. The 2012 CODATA evaluation Hi] combining e-p scattering and atomic transitions in hydrogen yields 
te = 0.8775(51) fm, 4% higher than the extremely precise muonic hydrogen result iH] which yields ve = 0.84087(39) fm, 
and inconsistent at the la level. The exact results from several of these extractions, along with our updated recommendation, 
are given in Table Uin Sec. |IV] 

Significant work has gone into attempting to understand the proton radius puzzle, examining corrections to atomic transitions 
in hydrogen and muonic hydrogen, as well as proposing explanations in terms of new physics. For detailed reviews, see Ref. M 
Till . In this work, we examine recent extractions of the proton radius from e-p scattering data and suggest a combined result 
and uncertainty that accounts for our current understanding of the systematic uncertainties and model dependence in such 
extractions. This will allow for a more reliable comparison of the scattering results, the atomic hydrogen measurements, and the 
muonic hydrogen results. 


II. EXTRACTION OF THE RADIUS FROM E-P SCATTERING 


Elastic e-p scattering depends on the proton charge and magnetic form factors Ge(Q 2 ) and Gm(Q 2 ) which, in the non- 
relativistic approximation, are the Fourier transforms of the charge and magnetization densities. In this approximation, the RMS 
radius of the charge density is related to the slope of the charge form factor at Q 2 = 0. Relativistic effects yield model-dependent 
corrections which invalidate this simple connection between the charge density and form factors. For the extractions examined 
here, the conventional approach was adopted, where the radius is defined in terms of the low-Q 2 expansion of the form factors 
as 


Ge(Q 2 ) — 1 — 'o.^eQ 2 + > (1) 

o 

allowing measurements of the form factors at low Q 2 to be used to extract the radius. In the one-photon exchange approximation, 
the spin-independent cross section depends on G 2 E and G 2 M , with the angular dependence at fixed Q 2 allowing for separation of 
Ge and Gm- In cases where one form factor dominates the cross section, spin-dependent cross section data, which are sensitive 
only to the ratio Ge/Gm, provide improved extraction of the form factor with smaller contributions to the spin-independent 
scattering lfl2l - fl4ll 

Most elastic e-p experiments provide a set of cross sections and uncertainties along with an estimate of the overall normal- 
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ization uncertainty which applies to the entire data set. A global analysis will often allow each of these normalization factors to 
vary, with a contribution to the fit’s total x 2 associated with the size of the normalization adjustment. The recent Mainz data 0 
has 31 normalization factors which are are completely unconstrained and must be determined in the fit, with the extrapolation to 
the known value Ge( 0) = 1 providing the absolute normalization. 

Within the world data set, excluding the recent Mainz measurement 0 , the individual experiments show consistency within 
the quoted errors SHU. However, a comparison of the Mainz form factors with world data clearly indicates that the Mainz data 
are not entirely compatible with the previous measurements, complicating the task of presenting a combined result. In addition, 
the way in which the Mainz data set quotes uncertainties is different from essentially all previous measurements. Simply treating 
the quoted uncertainties on the Mainz and world cross section measurements on equal footing would strongly bias the result 
in favor of the Mainz data, as discussed in more detail below. Because of this, we will examine the radius extraction from the 
Mainz data set and other world data separately, and then compare the results to examine the question of consistency. 


A. Mainz data set 

First, we note that because of the large number of data points in the Mainz data set 0 , coupled with the lack of measured 
absolute normalization of the cross sections, the A1 collaboration presented the data as several data subgroups, each with its 
own normalization, and separately estimated uncorrelated and correlated systematic uncertainties. The uncorrelated systematic 
uncertainties account for possible errors that are independent for every cross section and thus are treated identically to the 
statistical uncertainties. The correlated systematic uncertainties account for potential corrections which have a correlated effect 
across many measurements (but excluding pure normalization uncertainties). They chose to model the correlated systematics by 
applying corrections which depend linearly on the scattering angle, and which typically vary by 0.2% between the largest and 
smallest angle of each data subgroup. 

The separation of the uncertainties into uncorrelated, correlated, and multiple normalization uncertainties makes the extraction 
of the radius from the data rather complex. The data as provided in Ref. 0] already applies normalization factors determined in 
the global fit to the Mainz data, and no information on the original normalization is available. Therefore, any analysis of these 
data requires that the normalization factors of each subset be allowed to vary as part of the fit, as well as separate treatment of 
the uncorrelated and correlated systematic uncertainties. Using fixed normalization factors rather than allowing them to float 
as part of the fit is not consistent with the way the uncertainties are estimated by the experiment and can artificially decrease 
uncertainties by a factor of 5 or more 01. _ 

Several values of the charge radius have been extracted based on the Mainz data 0U1 EMI, and these analyses yield a 
significant range of results for te depending on the treatment of systematic uncertainties and normalization factors, choice of 
two-photon exchange corrections, and choice of fitting function used. However, most of these do not fully treat the uncertainties 
at the level required to obtain a reliable estimate of the uncertainty. The only complete treatments published so far are those of 
refs. 0E1, though another complete analysis which also includes a detailed examination of the systematic uncertainties has 
recently been completed Jlfj]. Ref. fH yields very different results for model-independent fits using the 2 -expansion l20fl and 
a dispersion relation analysis, and while their dispersion relation approach yields a radius that is consistent with the muonic 
hydrogen radius, it provides a much poorer fit to the Mainz data. Because of this, we will consider only the radius extractions 
shown in Ref. 0 in our examination of the Mainz data, although we will examine the assumptions and error treatment used in 
these extractions. 

As noted above, the Mainz breakdown of the uncertainties makes it difficult to perform a combined analysis of Mainz cross 
sections and other world data. In such fits, the relative weighting of the different measurements is determined by the uncertainties 
applied to each data point. For the Mainz data, this is only the statistical and uncorrelated systematic uncertainty, which repre¬ 
sents a small fraction of the total uncertainty, while other experiments include everything except the normalization uncertainty 
in these errors, including those associated with radiative corrections and two-photon exchange, neglected in Ref. 0 . Thus, in a 
naive global cross section analysis, the Mainz data will receive far more weight than it should compared to all other experiments. 


B. World data 

Data from several scattering experiments are available, e.g. see Refs. liTEi - jTHl [21~1 1221 , and there have been many parameteri- 
zations made of the proton form factors based on these data. In several cases, the goal was J;o have a global parameterization of 
the full data set, up to Q 2 m 30 GeV 2 for Gm- Commonly used global parameterizations I12ll I23l - l25il did not attempt to ensure 
that the fit function was flexible enough to provide a reliable extraction of the radius in the presence of a large body of higher-Q 2 
data. One exception is Ref. |[26j which performed a low Q 2 fit to constrain the slope, and then performed the global fit with the 
slope parameters fixed. Because we are not interested in the form factor at high Q 2 , we will examine only fits with the explicit 
aim of extracting the proton radius from low-Q 2 data. 
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In determining the radius from world data, we rely on the results from two of the more recent extractions SH which focus on 
the low-Q 2 region, examine the model dependence of the extraction, and use a relatively complete data set including polarization 
data. There are several other analyses which could be included but which we do not use here. Single-parameter fits are excluded 
as the Q 2 range necessary to provide a useful limit on the radius goes beyond the region where the single parameter fits are 
sufficiently precise to represent the data B, i27ll . Some fits |2j,;6L l28[f29 ll are excluded because they are essentially earlier versions 
of the fits which we do include, and as such do not add much new information. Others are excluded because they include only 
a limited set of the low Q 2 data, they_apply two-photon exchange corrections based on phenomenological extractions which 
are not well constrained at low Q 2 [30], or because they fit to extracted form factors rather than the original cross section 
measurements lt20l,[3lll which generally means neglecting normalization uncertainties on the different data sets. 


C. Two-photon exchange corrections 

One critical issue is the treatment of two-photon exchange corrections [32, 33] in the extraction of the charge radius ||28l 
134L [35]. The Mainz extraction 0 applies the correction of Ref |36], which is valid (up to corrections of order (Za) 2 ) only 
for scattering from a point nucleus. For the proton, the higher order corrections are negligible and in the limit Q 2 —> 0, the 
scattering is not sensitive to the internal structure, and thus the point-proton approximation is appropriate. At non-zero Q 2 
values, the correction changes in both the Coulomb distortion correction for finite-size nuclei [ 371 and the extension of hard TPE 
corrections to a proton with internal degrees of freedom, yielding a Q 2 -dependence II32L l33l 135L 38] which is clearly relevant in 
attempting to determine the slope of Ge(Q 2 ) as Q 2 —>■ 0. 

Note that the correction at Q 2 = 0 has the opposite sign and angular dependence compared to what is necessary to explain 
the discrepancy between Rosenbluth and polarization form factor measurements at larger Q 2 values Il2i I13L 1 1 51 122L 132l I39II . 
Calculations at finite Q 2 which account for the structure of the proton show this change of sign in going from Q 2 = 0 to larger 
Q 2 values [381. 14fj|- l42ll . but there is a significant model dependence associated with the modeling of the intermediate hadronic 
state. Recent comparisons of electron-proton and positron-proton scattering confirm a change of si gn by G 2 = 1 GeV 2 E 43l[44 1 
and are in reasonably good agreement with the TPE calculations performed in a hadronic basis I 38u40l - l42ll . In addition, all TPE 
calculations which are expected to be valid at low Q 2 values are in very good agreement below Q 2 = 0.2 GeV 2 [35]. Given 
the consistency of the calculations and the experimental demonstration that there is a significant Q 2 dependence in the TPE 
effects at low Q 2 , we will use only extractions which include hadronic TPE corrections, and include an estimate for the model 
dependence of these corrections. 


III. EXAMINATION OF PREVIOUS CHARGE RADIUS EXTRACTIONS 


In the following sections, we examine the existing extractions of the proton radius presented above and provide recommen¬ 
dations for updated radius values from Mainz and world data based on the inclusion of additional uncertainties and corrections. 
We also provide a simple combined result as well as a suggested combined result with additional uncertainty intended to account 
for the observe discrepancies between the Mainz and world data sets. 


A. Radius extraction from Mainz data 

We chose to take the extracted values of charge and magnetic radii from the analysis of Ref. 0 which includes the hadronic 
TPE corrections of Ref. fJsil : 

TE — 0.875(5) s tat (4)syst (2)model(^)group = 0.799(13) s tat (9)syst (5)model(3)group ; (2) 

where the breakdown of the uncertainty is detailed in Ref. 0 . Note that it is not clear in 0 ] how these uncertainties should be 
combined. While most authors combine these contributions in quadrature, yielding SrE = 0.008 fm, we note that all figures 
in Ref. 0 that show the combined uncertainties on the radii or form factors combine the errors linearly, corresponding to 
SrE = 0.016 fm. When asked about this, the response from the Mainz A1 collaboration was that the first three uncertainties 
should be combined in quadrature while the final contribution should be added linearly [45], yielding te = 0.875(12) fm, 
r M = 0.799(20) fm. 

As noted earlier, the quoted radius includes no uncertainty associated with the TPE correction. The impact of replacing the 
McKinley-Feshbach correction |36]] with the hadronic calculation [38] was small, decreasing te by 0.004 fm and increasing r m 
by 0.022 fm. We take half of this shift as an additional uncertainty to account for the model dependence of the TPE corrections. 

Refs, lfl^.1461 demonstrate that impact of correlated systematic errors on the radius can be significantly larger if the correction 
is assumed to have a different functional form different from that used in 0 (linear with 6) or if it was applied to different 
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subsets of the data (e.g. to each of the 34 data subgroups with different normalization factors rather than each of the 18 
spectrometer-beam energy combinations). Based on these tests, and concerns about other sources of correlated systematic 
uncertainties |5| 7f| (46] that are neglected in the extraction, we double the systematic uncertainty from on r E from 0.004 to 
0.008 fm, and on tm from 0.009 to 0.018 fm. 

With the additional TPE uncertainty, the increased correlated systematic applied to both charge and magnetic radii, and using 
the proposed procedure for combining the different sources of uncertainty, we obtain the following radii from the Mainz data: 
r E = 0.875(15) fm, r M = 0.799(28) fm. 


B. Radius extraction from world data 

For radii extracted from world data, excluding the recent Mainz cross sections, we take the extractions from Ref. H, r E = 
0.875(10) fm, and Ref. H, r E = 0.886(8) fm. Note that the result quoted in Ref. iH] was from a combined analysis of Mainz 
and world data, with additional systematic uncertainties applied to the Mainz data to account for the uncertainty in the subtraction 
of the target window contribution. The combined analysis yielded a radius that was larger by 0.001 fm than the result quoted 
here which excludes the Mainz data. 

These analyses include very similar data sets, with differences mainly associated with the choice of fitting function and the 
evaluation of the model dependence. We choose to take an unweighted average of the two results, r E = 0.881(9), and apply 
an additional uncertainty of 0.006 fm, equal to half of the difference between the two results, as an estimate of possible model 
dependence beyond what is included in the analyses. For the magnetic radius, the value of J 3 ]] is taken, as Ref. jj] did not provide 
a value of ?’m while the previous work jfj did not include the more recent polarization measurements which help constrain Gm 
at low Q 2 . This approach described above yields a final world data radii of r E = 0.881(11) fm and rju = 0.867(20). 


C. Combined e-p scattering radius 

Simply combining the charge radii from Mainz and world data yields r E = 0.879(9) fm. However, while the results for the 
charge radius from the Mainz and world data sets are consistent, there are significant tensions between the Mainz and world 
data, e.g. Figs. lOand 19ofRef. 0], In addition, the Mainz and world updated magnetic radii differ by three standard deviations 
as quoted in Refs. HU, and two standard deviations in our updated result (due to a combination of a slight reduction in the 
inconsistency and an increase in the quoted uncertainties). Given the tension between the results, it is not clear how one should 
treat the uncertainty on the combined result, and we are forced to make a relatively arbitrary decision on how to do so. 

We choose take the discrepancy in tm and use this to estimate the potential systematic errors in r E , based on the relative 
sensitivity of r E and r« to changes in the data. Examining the variation of the radii when applying corrections (e.g. TPE 
contributions) or varying fit functions, we find that the magnetic radius has a sensitivity that is roughly 10 times larger than 
that of the charge radius. We thus take 10% of the inconsistency as an additional systematic on the combined charge 
radius, yielding a final combined result for the charge radius of r E = 0.879(11) fm. For the magnetic radius, we simply 
take the weighted average and apply half of the difference between the Mainz and world results (0.034 fm) as an additional 
uncertainty. For the charge radius, this additional uncertainty yields a 20% increase in the total uncertainty, while it is the 
dominant uncertainty in the extraction of the magnetic radius. 


IV. CONCLUSIONS 

We provide recommended radii and uncertainties based on electron scattering measurements, taking into account the tension 
between different data sets. Note that the previous CODATA evaluation fH] included the extraction of Ref. d, which did not 
include the hadronic TPE contribution, and an earlier evaluation of world data 0 which had a smaller data set and larger 
uncertainty than the extractions of d 5i]. While we include additional data compared to the previous evaluation, our final 
uncertainty is noticeably larger as additional uncertainties have been applied to account for more detailed evaluation of some 
uncertainties and to account for tension between different extractions. 

Because of this tension, we are forced to make a somewhat arbitrary choice in evaluating the combined uncertainties, and have 
attempted to be somewhat conservative in the uncertainties we choose. Further examination of possible systematic uncertainties 
is presented in Refs. [5| Q3, but at this point the origin of the discrepancy between different data sets is not yet clear. 
We have laid out our assumptions and reasoning for the additional uncertainties included but without an understanding of 
the source of this discrepancy, it is not possible to come up with a more rigorous way of evaluating the combined result and 
uncertainty. We consider the approach presented above to be a reasonable combined average of all e-p scattering data, and even 
with this significantly more conservative approach, we still find a 3.5cr discrepancy between e-p scattering and muonic hydrogen 
measurements and a 5.1a discrepancy when combining the electron scattering and atomic hydrogen transitions. 
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Source 

te 

[fm] 

tm 

[fm] 

Published results 



pa [9] 

0.8409(4) 

0.870(60) 

eH [8J 

0.8758(77) 

- 

Mainz A1 [7,45] 

0.8790(110) 

0.777(19) 

Zhan [3] 

0.8750(100) 

0.867(20) 

Sick [5,6] 

0.8870(80) 

0.855(35) 

CODATA 12 average [8] 

0.8775(51) 

- 

New updates 



Mainz updated 

0.8750(150) 

0.799(28) 

world updated 

0.8810(110) 

0.867(20) 

naive global average 

0.8790(90) 

0.844(16) 

suggested global average 

0.8790(110) 

0.844(38) 


TABLE I: Charge and magnetic radii from various published extractions and from the averaging procedures described above. For ease of 
comparison, all charge (magnetic) radii are quoted to 4 (3) significant digits. 
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